Memory function on turbulent fluctuations in soft-mode turbulence 
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Modal relaxation dynamics has been observed experimentally to clarify statistical-physical properties of soft- 
mode turbulence, the spatiotemporal chaos observed in homeotropically aligned nematic liquid crystals. We 
found a dual structure, dynamical crossover associated with violation of time-reversal invariance, the corre- 
sponding time scales satisfying a dynamical scaling law. To specify the origin of the dual structure, the memory 
function due to non-thermal fluctuations has been defined by a projection-operator method and obtained numer- 
ically using experimental results. The results of the memory function suggest that the non-thermal fluctuations 
can be divided into Markov and non-Markov contributions, the latter is called the turbulent fluctuation (TF). 
Consequently, the relaxation dynamics is separated into three characteristic stages: bare-friction, early, and late 
stages. If the dissipation due to TFs dominates over that of the Markov contribution, the bare-friction stage 
contracts; the early and late stages then configure the dual structure. The memory effect due to TFs results in the 
time-reversible relaxation at the early stage, and the disappearance of the memory by turbulent mixing leads to a 
simple exponential relaxation at the late stage. Furthermore, the memory effect due to TFs is shown to originate 
from characteristic spatial coherency called the patch structure. 

PACS numbers: 61.30.-v, 05.45.-a, 47.54.De, 05.40.-a 



I. INTRODUCTION 

Weak nonlinearity can generate spatial and temporal disor- 
ders in systems where the number of effective degrees of free- 
dom increases with increasing system size. A phenomenon 
triggered by weak nonlinearity in such high-dimensional sys- 
tems is called spatiotemporal chaos, in contrast with chaos 
where unpredictable behavior emerges from few degrees of 
freedom in a deterministic way. Theoretical work on the spa- 
tiotemporal chaos (Refs. 01 Q] and references therein) has 
covered topics such as the complex Ginzburg-Landau equa- 
tion, the Kuramoto-Sivashinsky turbulence (KST), the Niko- 
laevskii turbulence, and coupled map lattices. An outstanding 
feature of chaos and turbulence is its dual structure, where the 
dynamics is separated into an initial regime corresponding to 
deterministic short orbits and a final regime corresponding to 
stochastic long orbits. Mori and Okamura, for example, have 
theoretically studied 1D-KST and revealed the dual structure 
in turbulent mixing by numerical analysis yfl. 

Convection systems have been experimentally investigated 
to study nonlinearity. An example is electrohydrodynamic 
convection observed in the nematic liquid crystal and con- 
trolled by an ac voltage applied to system. One can easily 
achieve accessible characteristic length and time scales in ex- 
periments of electroconvection compared with the Rayleigh- 
Benard convection. Another advantage of investigating ne- 
matics is the ease in controlling the anisotropy to study sym- 
metry properties. There are two types of layer alignment 
in nematic liquid crystals; one is planar alignment in which 
the director aligns parallel to substrates (x-direction), and the 
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other is homeotropic alignment in which the director aligns 
perpendicular to substrates (z-direction). Rubbing along the 
x-direction of a substrate's surface (x—y plane) produces pla- 
nar systems, and intrinsically breaks the rotational symmetry. 
In homeotropic systems, by contrast, the rotational symmetry 
in the x—y plane remains. With a sufficiently strong applied ac 
voltage, magnitude V and fixed frequency /, the Freedericksz 
transition, occurring at a certain threshold voltage Vf, spon- 
taneously breaks the rotational symmetry J4J]; the transition 
is accompanied by the excitation of the Nambu-Goldstone 
modes Q5j— LZD ■ By further increasing V, electrohydrodynamic 
convection occurs at V c . The nonlinear coupling between the 
convective and Nambu-Goldstone modes generates a pattern 
that is both spatially and temporally disordered. The experi- 
mentally observed phenomenon, called soft-mode turbulence 
(SMT), is an example of spatiotemporal chaos ifsl— T3~oh . 

In our preceding study, we have observed SMT relaxation 
by measuring a temporal autocorrelation function and re- 
ported that the relaxation is well fitted by a compressed ex- 
ponential function ifTTll . Because the compressed exponen- 
tial is employed to describe the dynamics of jammed systems 
JI2I1 . we remarked on the similarity between SMT and glass 
forming liquids. The nonlinearity in the dynamics of the lat- 
ter originates from dynamic coherency in some regions. In 
fact, spatiotemporal fluctuating cooperative regions have been 
observed as dynamical heterogeneities near the glass transi- 
tion point [13-16]. A characteristic length for these coop- 
erative regions increases as the glass transition point is ap- 
proached. In SMT, on the other hand, patch domains exist in 
which convective rolls align in a unique orientation lfl7l - fl9ll . 
The characteristic size £ of a patch domain is several times 
longer than the typical size of convective rolls and decreases 
with the distance from SMT onset; £ ~ e^ 1 / 2 with control 
parameter e. We have therefore concluded that the SMT patch 
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domains behave like the cooperative regions in glass forming 
liquids and the coherent motion in the domains generates non- 
exponential relaxation. To study SMT dynamics in detail, the 
temporal correlations of each wave number, modal autocorre- 
lation functions, are suitable. Our previous study focused on 
the net autocorrelation function consisting of the entire wave- 
number information. Here, we investigate the modal relax- 
ation dynamics to specify statistical physics of SMT. 



H. EXPERIMENT 

We study a 2D pattern dynamics of SMT observed in a 
homeotropic alignment of nematic liquid crystals. This study 
follows a standard setup [H [2(J, 12111 . The space between 
two parallel glass plates, spaced 27 fj,m apart, was filled 
with the nematic liquid crystal, iV-(4-Methoxybenzilidene)- 
4-buthylaniline (MBBA). The plate surfaces were coated with 
transparent electrodes, made of indium tin oxide (ITO) with a 
circular cross-section of radius 13 mm. To obtain homeotropic 
alignment, the surfaces was covered by a surfactant, N, N- 
dimethyl-iV-octadecyl-3-aminopropyl-trimethoxysilyl chlo- 
ride 50% (DMOAP). The values of the dielectric constant and 
electric conductivity parallel to the director were en = 6.25 
and = 1.17 x 1Q - ' fi -1 m — respectively. Denoting the 
dielectric constant perpendicular to the director by ej_, the di- 
electric constant anisotropy e a = er|| — ej_ is found to be neg- 
ative. An ac voltage V(t) — \/2V cos(2tt ft) was applied to 
the sample. For a control parameter, we employ a normal- 
ized voltage e = (V/V c ) 2 — 1, where V c denotes the threshold 
voltage for electroconvection having value 7.78 ± 0.05 V. We 
show results fore = 0.025, 0.050, 0.075, 0.10, 0.20, 0.30, and 
0.40, where e has a margin of error of ± 0.013. Another con- 
trol parameter was the frequency / of the ac voltage. Two 
patterns of SMT arise; oblique roll in / < /l and normal roll 
in / > /l, where /l denotes the Lifshitz frequency We 
set the frequency / = 100 Hz well below /l. The temperature 
was regulated to 30.00 ± 0.05 °C. Before each sampling, we 
waited for 10 min at fixed voltage V w and then a further 10 min 
at the set V to get a desired e value, where Vf < V w < V c . 
The waiting time is sufficiently long for systems to achieve 
steady state. 

The electroconvection pattern was observed by a micro- 
scope (ECLIPSE E600POL, Nikon Corporation, Tokyo) and 
was captured by a high-speed camera (HAS220, DITECT 
Co. Ltd., Tokyo) that can successively take 4200 frames. The 
numbers of frames per second were 5 for e < 0.1 and 10 oth- 
erwise. A typical two-dimensional image is shown in Fig. Q] 
(a). The measurement area was 830 x 830 fim 2 (450 x 450 pix- 
els). The transmitted light intensity I(x, t) at each pixel was 
digitized into 8-bit (i.e., 256-level) information, where a series 
of pattern analyzing was processed according to Ref. [ 22] . 

With angle brackets signifying the long-time average in the 
steady state, the temporal correlation of two functions is de- 
fined as 
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FIG. 1. Static information of our experiment at e = 0.1. (a) A 
typical snapshot of SMT with the white scale bar indicating 100 /jm. 
Each bright line indicates upward flow, (b) Magnitude tints of the 
spatial power spectrum Pk in the Fourier space, where shading cor- 
responds to the Pfc value, (c) Plot of the power spectrum P k as a 
function of the radial wave number k = |fe|. The peak, marked by an 
arrow, corresponds to the fundamental period in a convective rolling. 



We employ the Fourier transform of the fluctuation 

AI(x,t) = I(x,t) — (I(x,t)} as the gross variable Wfe(t); 



u k {t) := J dxAI{x 7 t)e ikx , 



(2) 



where i = and the integral range is over the entire sys- 
tem domain. In isotropic systems, it is sufficient to study u k (t) 
1I23I1 . where k denotes the radial wave number; k = \k\. We 
focus on the normalized modal time-correlation function of k; 



U k (r) := (u k (t + r)u* k (t))P k -\ 



(3) 



(f(t + r)g(t))= lim -L 

T— ¥00 11 



dtf(t + T)g(t). (1) 



where the asterisk denotes the complex conjugate operation 
and P k = (\u k (t)\ 2 ') denotes the spatial power spectrum. 

Note that U k (r) is a real number due to translational sym- 
metry and isotropy. The spatial power spectra as functions 
of the wave number k and k = \k\ are illustrated in Fig. Q] 
(b) and (c), respectively. SMT isotropy is clearly reflected in 
Fig.rjJ(b). Although experimental parameters in Ref. lf23ll are 
slightly different from ours, the P k profiles agree. A clear 
peak exists at fc pea k — 0.321 ^m -1 for e = 0.10 and its cor- 
responding length Ap ea k = 27r/fc pe ak is 19.6 /im that is a half 
of Ao, where Ao = Ao(e) denotes the length of a pair of elec- 
troconvections. The peak wavenumbers /c pea k for other e val- 
ues are similar. The wave number k is normalized by Ao as 
k = kAo/2Tt. 
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III. RESULT AND DISCUSSION 



decay form; 



A. Dual relaxation 



U k (r) cx 1 



/ (a) 



(4) 



The modal autocorrelation function (t) obtained exper- 
imentally is shown in Fig. In the short-time regime, the 
relaxation is not described by a simple exponential function. 
Instead, we found dual relaxation as evident in Fig. [3] In the 
early stage [24], the correlation function follows an algebraic 




In contrast, the relaxation in the late stage is well described by 
exponential decay; 



U k (r) cx exp 



(5) 



Figure@]presents the characteristic time scales at each stage 
for the dual structure. The time scales are almost constant at 
small wave numbers. Note that the results at small k may be 
affected by limitations in sample averaging. At large wave 
numbers, the time scales of the dual relaxation satisfy dy- 
namic scaling laws; 



oc k 



t% cx k 



(6) 



with dynamical exponents z a ~ 1.0 and z e ~ 1.5 regardless 
of e. Although Mori and Okamura have reported that the dy- 
namical scaling exponents obtained theoretically in 1D-KST 
are z. d = 1 and z e = 2 Jl, we cannot directly compare our 
results with their theoretical work because the dynamical ex- 
ponent depends on spatial dimensionality. Indeed, dynamic 
scaling for several turbulence models depends on the dimen- 
sionality II 25142711 . 

To elucidate the dual structure of SMT theoretically, we 
have taken account of an experimental fact that the spatiotem- 
poral disorder accompanied by the patch structure can be rep- 
resented as non-thermal fluctuations lfl7[|28ll . The projection- 
operator method for turbulence rf29ll is one of the most suitable 
methods to manage such non-thermal fluctuations; its math- 
ematical procedure derives an evolution equation containing 
the memory effect that is equivalent to the autocorrelation of 
fluctuations (see Appendix|A]i. Thus, in the rest of the paper, 
analysis for the memory functions plays a central role to study 
a relationship between the dual structure and the non-thermal 
fluctuations. 




FIG. 2. (color online) Plot of the modal correlation functions at e = 
0.050 (top), 0.10 (middle), and 0.40 (bottom). The wave numbers 
k are 0.76 (red square), 1.0 (brown inverted triangle), 2.0 (black 
circle), 3.0 (green triangle), and 4.0 (blue diamond). 



FIG. 3. (color online) The normalized modal autocorrelation func- 
tion of e — 0.1 and k = 2.0, which is the same as shown in Fig. [2] 
The algebraic decay l|4j (brown solid line) well describes the dynam- 
ics in the early stage, but a simple exponential l[5} (green dashed line) 
is better at the late stage. 
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FIG. 4. (color online) Log-log plots of time scales T k (top) and 
r k (bottom) for several control parameters; e = 0.025 (pink), 0.10 
(blue), 0.20 (green), 0.30 (brown), and 0.40 (red). The dashed lines 
plot the power law with exponents z a = 1.0 (top panel) and z e — 1.5 
(bottom panel). 



B. Memory effect due to non-thermal fluctuations 

We assume the modal elements {u k } form a complete set 
of macroscopic variables in SMT 1 30] . According to the 
projection-operator formalism, the evolution equation in SMT 
is represented as 



dU k (r) 
dr 



dr'T' k (r-T')U k (r'), 



(7) 



where T' k (r) denotes the memory function that results from 
the non-thermal fluctuations. Here, translational symmetry re- 
duces the temporal correlation of the modal elements to 



(u h (t + T)u* k ,(t)) = 6 k , k > (u h (t + T)u* k (t)) 



(8) 



and the mechanical coefficients ui kk are zero by definition. 

Using experimental results for the modal correlation func- 
tion, one can numerically solve Eq. (0 to obtain the memory 
function. The memory function T^r) has a sharp peak at 
r = 0, as evident in Fig. [5] implying that the non-thermal 
fluctuations can be separated into Markov and non-Markov 
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FIG. 5. The memory function T' k (r) characterized by Eq. 10. This 
is the result of e — 0.1 and k — 2.0 as an example. 



contributions. Therefore, with 



T / k (T) = 2 1 ^S(r)+T k (r), 



(9) 



where 7[°' denotes the bare friction due to the Markov contri- 
bution of the non-thermal fluctuations and T k (t) the memory 
function due to the non-Markov contribution, the evolution 
equation (0 reduces to 



dr'r fe (r - T')U k {T'). (10) 



To emphasize the transport due to the turbulent-like dynamics, 
the non-Markov contribution to the non-thermal fluctuations 
is here called the turbulent fluctuations (TFs) in SMT. One 
can define a characteristic time scale for the memory function 
due to TFs as 



_(r) ._ 



1 



r fc (o) 



diT fc (T), 



(11) 



within which the memory effects due to TFs are alive. 

As seen in Fig. [6] the memory function Tk if) caused by TFs 
has non-negligible magnitudes and time scales over a range of 
wave numbers in the early stage; therefore, SMT dynamics is 
regarded as being non-Markovian due to TFs. The memory 
has a small peak around r ~ 2 s. We believe that this is 
due to a residual of the mechanical coefficient appearing in 
the projection-operator method, despite the coefficient being 
theoretically zero. 

Each equation of the dual relaxation is elucidated by the 
evolution equation (llOt . First, let us consider the dynamics at 
the late stage, 7V r. The frequency-dependent friction 
Tk, u caused by TFs is defined as the Fourier-Laplace trans- 
formation of the memory function; 



drr fc (r)e k 



(12) 




FIG. 6. (color online) Time-dependence of the memory function for 
e = 0.10 at several fixed wave numbers: k = 0.76 (red square), 1.0 
(brown inverted triangle), 2.0 (black circle), 3.0 (green triangle), and 
4.0 (blue diamond). 



At the late stage, friction can be regarded as static, 

tF := T k , u=0 ; (13) 
hence, the memory effect is approximated by a delta function 

r fe (r) ~ 2jPs(t). (14) 

The modal relaxation dynamics is thus obtained as a simple 
exponential function (O with the characteristic time 

^ = (7f +7 f)". d5) 

Next, at r <C , an approximate solution of Eq. (TTOb is 

(r) r r 

Mr) ^ 1 - ^frj 7 o dT ' (Vf ' + r') + 0(t 3 ) (16) 

with 7^ = 7fc /7I ■ Thus, if % <§; 1, the time range can 

(r) (r) (r) 

be separated into r <C JkTt and 7 fcT£ <C r <C 77 . In 

(r) 

the shorter time regime r <C 7fcT^ , called the bare-friction 

stage, modal relaxation decays linearly with slope 7^. In 
contrast, it reduces to time-reversible algebraic decay in the 

(T) (D 

longer time regime 7^.77 ^ T ^ T k ' i- e -' trie eal "ly stage 
of the dual structure with characteristic time 

4* = f¥T¥- (17) 

Therefore, a small ratio % of the friction coefficients is a nec- 
essary requirement for the appearance of dual relaxation in 
SMT. The characteristic stages are summarized in Table U 



FIG. 7. (color online) Wave-number dependences of the time scale 
r k for several values of e; e = 0.050 (purple), 0.20 (green), and 
0.40 (red). Each dotted curve is obtained by a peak-fitting method 
and the vertical arrow points where kr is. 

C. Characteristic time scales 

A characteristic time scale for the modal autocorrelation 
function can be defined by 

/>oo 

rf ) := / dr^(r). (18) 
Jo 

This is the relaxation time due to non-thermal fluctuations. 
In contrast, the time scale 1 /7I denotes the relaxation time 
due to Markov fluctuations only. From the evolution equation 
([Tol l, is exactly represented by the friction coefficients, 

rf = (7f +7F)" 1 . d9) 

Therefore, is nearly equal to according to Eq. H5h 
its accuracy has been confirmed within the margins of numeri- 
cal error, indicating that the approximation ( fl~4l > is reasonable. 

The characteristic time scale 77. , plotted as a function of 
normalized wave number in Fig. [7] has a weak peak in k < 2 

(V) 

for each e. Let fcr denote the wave number at which 77 ' 
has a peak, where fcr was determined by a method for multi- 
peak fitting so that roughness of curves was smoothened. In 

expectation that the coherency in patch domains leads to the 

(r) 

non-Markov memory effect, the time scale 77 should in- 
clude features relating to the patch domains. The length scale 
Ap = 27r/fcr is several times larger than the diameter of an 
electroconvective roll. In addition, the power law Ar oc e -1 / 2 
is quantitatively reasonable as indicated in Fig. [8] Therefore, 
it follows that the dual relaxation caused by TFs originates 
from the patch structure. The time scale rapidly decreases at 
larger length scales fc < fcr, indicating inter-patch dynamics 
does not affect the memory effect due to TFs. Some work 
does support the relationship between the patch structure and 
the dual structure B17ll28ll31ll . where the dynamics of a tagged 
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TABLE I. Characteristic stages of modal relaxation dynamics in soft-mode turbulence. 



stage 


time range 


decay form 


Bare-friction stage 


T « %ri n 


Linear decay caused by the Markov fluctuations 


Early stage 


%rP « r « rP 


Time-reversible decay that originates from memory due to TFs [Eq. {4}] 


Late stage 


rP « r 


Simple exponential decay after turbulent mixing [Eq. lf5](] 



0.04 



0.02 



0.00 




0.0 0.1 



FIG. 8. e-dependence of fcr at which r k has a local maximum. The 
gray dashed line marks the line fit, fcp = co + cie with co = 0.0045 
and ci = 0.052. The exponent of the power law is suggested from 



the SMT patch structure: £ ~ e 



-1/2 



particle in the SMT disorder can be divided into two types of 
mode: one dominated by convective rolling within patch do- 
mains (i.e., intra-patch dynamics) and the other dominated by 
transfers with slow patch movements (i.e., inter-patch dynam- 
ics). 

Figure [9] shows corresponding dispersion relations. The 
friction 7^ is one order larger than the bare friction near 
kr, but converges at large wave numbers. A characteristic fea- 
ture is a dip appearing near k ~ 1 in but not in 7^ . The 



bare friction 7^ and the static friction 7^ caused by TFs 
are explicitly related to the average rate in entropy production 
Sk 029I1 . Over macroscopic time scale, it is analytically repre- 
sented as 



S k 



(7i 0) + 7f)- fc B/, 



(u) 



(20) 



where fc B is the Boltzmann constant. The dip might reflect a 
law of minimum entropy production rate in electroconvection. 



The time scales, t£ , rjf\ and %t^', are compared in 

Fig.[l0]for several e values. The correlation time is much 

longer than at large length scales, where the relaxation 
dynamics is approximately represented by a simple exponen- 



(n 



tial. The ratio ri /ri approaches unity with increasing k, 
signifying that the memory effect persists for relatively long 
times at short length scales. Meanwhile, 7J.7V ; approaches 

at large k, implying that the duration of the early stage 
is shorter at small length scales. Therefore, dual relaxation 




FIG. 9. (color online) Dispersion relations for several e values; e = 
0.050 (purple), 0.20 (green), and 0.40 (red). The top panel gives the 



bare friction coefficient 7} 
caused by TFs. 



(0) 



and the bottom panel gives friction 7^ 



clearly appears in SMT at intermediate length scales, where 



Tu is not too short compared with t^' and the friction co- 
efficient due to TFs is much larger than the bare friction, i.e., 

7fc < 1- 

Mori and Okamura have proposed the affinity hypothesis 
between the correlation and memory function [ 3], in which the 
correlation and memory functions have the same form of de- 
cay with different magnitudes and characteristic time scales. 
The main prediction from the hypothesis is that, when ^ 

is almost the same as r k , decay of the correlation features a 
long-time tail with exponent —3/2 over all spatial dimensions. 
Those time scales for SMT are almost the same for k > 3.0, 
as represented in Fig.[l0] Our modal correlation functions in 
the long-time region however have a large margin of error; we 
cannot evaluate whether such long-time tails exist or not. 
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FIG. 10. Plot of the characteristic times r^ U ' (dashed line), 
(dark solid line), and 7^17, (light solid line) for e = 0.050 (top), 
0.10 (middle), and 0.40 (bottom), respectively. 



IV. CONCLUDING REMARKS 

We experimentally observed the modal autocorrelation 
function U k (t) of pattern dynamics in SMT. The modal re- 
laxation dynamics featured a dual relaxation (Fig.0; the cor- 
relation is well described by the time-reversible function at 
the early stage and by a simple exponential relaxation at the 
late stage. The corresponding time scales and rj^' obey 
dynamic scaling laws for the intra-patch scales (Fig.|4|i. 

With an expectation that the patch structure generates mem- 
ory effects due to the non-thermal fluctuations, we employed 



the projection-operator formalism in which the memory func- 
tion naturally arises as the autocorrelation function of the non- 
thermal fluctuations. Indeed, we have verified the expecta- 
tion; the peak wave number fcr (Fig- ID of the time scale t% 
(Fig- 13 suggests that the dual relaxation is caused from intra- 
patch dynamics. 

Solving the evolution equation [Eq. (0] derived in 
the projection-operator formalism, we revealed two well- 
separated contributions for the non-thermal fluctuations in 
Fig. [5] rapidly-varying Markovian fluctuations and non- 
Markov TFs. The former relates to the bare friction 7^.° and 
the latter to the memory function T k (t) [Eq. ©]. SMT dy- 
namics were shown analytically to separate into three stages 

where the 



(r) 



(Table JJ: the bare-friction stage, r <C -y^r^ 

(r) 

relaxation is linear decay; the early stage, Jk T k <?C r <C 

, where the relaxation is time-reversible algebraic decay 

[Eq. (0J]; and the late stage, *C t, where the relaxation is 
simple exponential decay [Eq. ||5)]. Further, by the variables 



defined in the theory (i.e. 



(0) 
Ik > 



JfFK and ), we have an- 

(») 



r k 



r, ' [Eq. ©]. 



alytically specified time scales obtained in experiments: t, 

[Eq. O], 4 e) [Eq. 03}], and 

Comparison between the characteristic time scales (Fig.fTUb 
clarified the nature of the dual structure in SMT. It clearly 
appears at intermediate length scales where the characteristic 
time r^ u ' of the modal relaxation is not too short compared 
with that r k ' of the memory function due to TFs. Also, the 



(T) 

friction coefficient 7^ due to TFs is much larger than the 

bare friction "f^ ■ 

The physical origin of the Markovian non-thermal fluctua- 
tions is still an open question. The dip near k ~ 1 in Fig. [9] 
suggests that the Markov contribution contains spatial features 
relating to electroconvection. That is, the bare friction con- 
sists of not only thermal fluctuations but also rapid variation 
in non-thermal ones. Note that the Markovian non-thermal 
fluctuations seem not to affect the SMT pattern dynamics, but 
acts as a trigger for the onset of SMT. If non-thermal Markov 
fluctuations played a role in pattern dynamics, the modal cor- 
relation would break time-reversal symmetry even at the early 
stage; nevertheless, early-stage relaxation is invariant under 
time reversal [Eq. (@)]. 

Besides SMT, several types of spatiotemporal chaos appear 
in the electroconvection of the liquid crystal systems. For ex- 
ample, instead of the homeotropic anchoring, one can apply 
the memory-function procedure to the planner anchoring as 
the boundary condition in which the rotational symmetry is 
not alive in the x-y plane. Our approach proposed in this pa- 
per can pave the way to understand them universally from a 
viewpoint of the non-thermal fluctuations relating to the pat- 
tern. This will be discussed elsewhere. 
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Appendix A: Projection-Operator Formalism 

Physically, a macroscopic system can be divided into 
slowly varying behavior described by a set of macroscopic 
variables {Ai} and rapidly-varying terms IT321 l33ll . The 
projection-operator formalism mathematically allows us to 
separate these two dynamics. 

Using the linear projection operator, Mori mathematically 
derived a generalized linear Langevin equation containing 
a memory function [34]. The memory function is repre- 
sented as temporal correlation of the fluctuations r.;(t), where 
(rj(t + r)Aj(t)) = 0. The generalized linear Langevin 
equation can be employed for not only equilibrium but also 
non-equilibrium systems. However, in the linear projection 
scheme, the fluctuation term Ti{t) is orthogonal only to the 
linear functions of {Ai] and can consist of not only the micro- 
scopic degrees of freedom but also fluctuations from nonlin- 
ear terms of {Ai}. Such non-thermal fluctuations are possibly 
relevant in nonequilibrium systems 11351 136TI . 

Zwanzig formulated the generalized nonlinear Langevin 
equation by the nonlinear projection operator l37ll38ll as 

A i {t) = Vi{t) + J i {t)+Rf\t), (AI) 

where Vi(t) denotes the streaming term including nonlinear 
reversible terms and Ji (t) an irreversible term. The fluctuation 



term J?> ' (t) satisfies 

(R ( ? ) (t + T)F(A(t)))=0 

for an arbitrary function J 7 (A) of a set of macroscopic vari- 
ables, {Ai}. The nonlinearity of {Ai} is included in Vi(t) and 
Ji(t), and one can extract the linear part of {Ai} from Vi(t) 
and Ji (t) using the linear projection-operator method. Mori 
and Fujisaka derived the evolution equation for the correlation 
function Uij(T) := {A t {t + r)A*{t)) as JUS 

(A2) 

where uj.ij denotes the mechanical coefficient and T' a (r) the 
memory function. Note that Eq. (IA2t . which relates the mem- 
ory r'(r) to the non-thermal fluctuations, is of the same form 
as the evolution equation derived in the linear projection- 
operator formalism; therefore, from a physical point of view, 
the obtained memory functions should be checked to see 
whether the memory effect originates from non-thermal fluc- 
tuations or not. If the characteristic time scale for the correla- 
tion of R\ f (t) is extremely short, then R\ f (t) satisfies; 

(R?\t + T)R?\t))^^5(T), (A3) 

where 7}? denotes the bare friction coefficient. In addition, 
if the fluctuations extracted from the nonlinear terms do not 
corr elate with R^\t), the memory function due to the non- 
thermal fluctuations can be divided into two different terms, 

T' ij (r) = 2 7 ifS(r)+T ij (r), (A4) 

where F^t) is regarded as the temporal correlation of the 
non-Markov contribution of the non-thermal fluctuations. 
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